Impact of layer defects in ferroelectric thin films 
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Abstract 

Based on a modified Ising model in a transverse field we demonstrate that defect layers in 
ferroelectric thin films, such as layers with impurities, vacancies or dislocations, are able to in- 
duce a strong increase or decrease of the polarization depending on the variation of the exchange 
interaction within the defect layers. A Green's function technique enables us to calculate the polar- 
ization, the excitation energy and the critical temperature of the material with structural defects. 
Numerically we find the polarization as function of temperature, film thickness and the interaction 
strengths between the layers. The theoretical results are in reasonable accordance to experimental 
datas of different ferroelectric thin films. 
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I. INTRODUCTION 



Defects in crystals influence significantly the physical properties of almost all materials. 
In the last decades defect engineering has been developed as a part of modern semiconductor 
technology. There is also an increasing interest in studying defects and their related strain 
fields in ferroelectrics (FE), in particular in low dimensional systems [if. A broad variety of 
defects have been implemented and analyzed mainly in AB03-type FE, such as BaTiO^and 
Pb_fZr,Ti)0 3 (PZT) [LAS, ^, bismuth laver-structured FE, such as Bi 4 Ti 3 0i 2 (BTO) flS 
Q3aswellasSrBi 2 Ta 2 9 (SBT)Q,H To control the polarization and piezoelectric 
properties, PZT is often modified by various elements of lower valency (K, Mn, Fe) and 
higher talent (La, Nb, Ta) cations. The modified PZT are generally classified into the 
two categories named as "hard" and "soft" ones Q|. The substitution of lower- valent 
cations, compared to the constituent ions, leads to oxide vacancies which induces a soft FE 
behavior, whereas higher-valent substitution gives rise to cation vacancies offering a hard 
FE behavior. Contrary to hard FE, the soft FE are characterized by a lower coercitive field 
E c , lower remanent polarization a r and higher hysteresis losses. The different FE properties 
are generally originated by the occurrence of defects. In materials with confined dimensions, 
the influence of defects can be even more pronounced due to the enhanced relative volume 

n 

of the "defective" regions |4|. 

FE thin films have been extensively studied due to their ability for applications in non- 
volatile FE random access memories (FeRAMs) Q]. For this application the thin film 
ferroelectrics should reveal large remanent polarizations, low coercitive fields, and fatigue- 
free properties. PZT and SBT have been widely investigated as appropriate material in 
Fe-RAMs. Otherwise the poor fatigue characteristic and the small remanent polarization, 
are viewed as the major barrier for their applications. 

BTO is a promising alternative to PZT, because of its large spontaneous polarization 
along the a-axis. Recently, it was reported that the substitution of small amounts of im- 
purities, such as Sm and Nd for Bi and V, W and Nb for Ti in the pseudoperowskite 
[(Bi 2 Ti 3 Oio) 2- ] layers of BTO, is effective for lowering the leakage current and enhancing 
the remanent polarization (f| Q, ^ . Noguchi et al. jj reported the control of FE proper- 
ties in Sr-site-modified SBT bulk ceramics, in which the Sr-site was partially substituted 
by La atoms, showed a lower coercitive field as the original material; but, the remanent 



2 



polarization was also decreased in this material. On the other hand, several reports have 
revealed that the a r value is enhanced in Sr-deficient and Bi-excess SBT films because of 
the substitution of excess Bi atoms in the Sr site ^3|. However, E c also increases in 
these films. Since the ionic radius of the La ion is almost the same as that of the Bi ion, the 
electronegativity is also important for controlling the FE properties in the site-engineering 
technique. In fact, the electronegativity of a Bi ion is much higher than that of a La ion. A 
similar effect was obtained in Sr-deficient and Pr-substituted SBT films [lo| . 

PZT films are considered to most promising candidates for applications in microelec- 
tromechanical systems (MEMS) since they have large piezoelectric coefficients and elec- 
tromechanical coupling coefficients. As known from bulk ceramics, the piezoelectric and 
dielectric properties vary substantially with the composition (Zr to Ti ratio). While the 
bulk PZT ceramics have been well investigated, there is still a lack of understanding the 
piezoelectric and dielectric properties of PZT thin films, particularly as a function of compo- 
sition. With increasing Ti-content, larger remanent polarization and higher coercive voltage 
are observed []|. Impact of misfit dislocations on the polarization instability of epitaxial 
nanostructured FE perovskites (PZT nanoislands) is investigated by Chu et al. |4(. Their 
results suggest that misfit engineering is indispensable for obtaining nanostructured FE's 
with stable polarization. The properties of BaTiC>3 are modified by substituting of tetrava- 
lent ions such as Zr 4+ for Ti sites. This increases the dielectric constant and broadens 
the sharp temperature dependence of the dielectric constant around the Curie temperature. 
With the addition of strontium to PbTiC>3, the phase-transition temperature decreases with 
the increasing Sr concentration. 

Influence of defects on FE phase transitions is one of key-like topics of the modern studies 
of FE. Constant attention to this topic is connected with its role in the real FE problem 
solution from the point of views of the fundamental science as well as of the above mentioned 
applications. It is well known that localized spin excitations can arise in FE and magnetic 
systems with broken spatial translational symmetry, i.e. due to the presence of boundary 
surfaces, interfaces, impurities and other defects in the crystal. These localized modes ap- 
pearing in the above mentioned structures are experimentally observed using the neutron 
and Raman scattering [lgI . \y\ , or, more recently, far infrared absorption . The Green's 
function theories have been extensively employed to describe these modes in semi-infinite 
transverse Ising models ^| and FE thin films j^J. In addition to bulk spin waves there 
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may be localized modes associated with the impurity layer and with the surface j^j]. 
A transfer matrix formalism is applied to study a semi-infinite Heisenberg ferromagnetic 
superlattice with surface cell and impurity cell at arbitrary distance from the surface j^lj. 
The occurrence of localized modes associated with the impurity cell is demonstrated suffi- 
ciently In the dielectric continuum approximation, Zhang et al. f.22] investigated the effect 
of different kinds of coupling to the localized interface optical-phonon modes in two coupled 
semi-infinite iV-constituent superlattices with a structural defect layer. The influence of 
molecular impurity ions on FE phase transitions is studied by Vikhin and Maksimova using 

n 

the phenomenological Landau theory [23]. It is shown that molecular impurity ions with 
charge transfer-local vibration degrees of freedom which bi-linearly interact with FE order 
parameter can induce rather strong increasing of the critical temperature of FE phase tran- 
sition at reasonable concentrations of molecular impurities. It could be topical, for instance, 
for the KDP:Mn04 case. Using a combination of first-principles and effective-Hamiltonian 
approaches, Dieguez et al. j3| obtained a phase diagram in temperature and misfit strain 
in epitaxial BaTi Oa that is qualitatively different from that reported by Pertsev et al. (2^. 
In particular, in |2jJ is found a region of "r phase" at low temperature where Pertsev et al. 
have reported an " ac phase" . 

The most of the published papers consider the defects influence on the spin excitations. 
A first attempt to find out the influence of defects on the polarization of FE thin films are 



elucidated by Alpay et al. 



271 ] . A thermodynamic analysis has been carried out to 



investigate the role of dislocations in FE materials. Due to the coupling of the stress field of 
the dislocation and the polarization, there is expected a drastic variation of the polarization 
near to the dislocation. The compressive regions enhance the polarization and increase 
Curie temperature whereas tensile stresses decrease polarization with a commensurate drop 



in the Curie temperature 



26|, El- But actually the dependence of polarization and Curie 



temperature on defects in FE thin films is not so intensively studied theoretically. Therefore, 
it is the aim of the present paper to investigate the polarization of FE thin films with defect 
layers in more detail. Based on the Ising model in a transverse field and using the Green's 
function theory we calculate the spin-wave energies, the polarization and the phase transition 
temperature for a FE thin film with different structural defect layers. The numerical results 
are compared with those of bulk materials and with thin films without defect layers. 
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II. THE MODEL AND THE MATRIX GREEN'S FUNCTION 



Let us consider a three dimensional ferroelectric system on a simple cubic lattice composed 
of iV layers in z-direction. The layers are numbered by n = 1, . . . , N, where the layers n = 1 
and n = N represent the two surfaces of the system. The bulk is established by the remaining 
(N — 2) layers. The specific surface effects are included by additional coupling parameters 
between bulk and surface layers. In particular, we start with the Hamiltonian of the Ising 
model in a transverse field which includes both, bulk and surface properties. 

H = -\ £ - - n > £ s ? - n * £ c 1 ) 

ij ieb ies 

where S x and S z are components of spin-^ operators, and fit, and fl s represent transverse 
fields in the bulk and surface layers, and the sums are over the internal and surface lattice 
points, respectively. Jy is an exchange interaction between spins at nearest-neighbor sites i 
and j, and = J s between spins on the surface layer, otherwise it is Jf,. We assume that 
one or more of the layers can be defect, since Jd and Qd denote the exchange interaction 
and the transverse field of the defect layer. In the ordered phase we have the mean values 
< S x >^ and < S z >^ 0, and it is appropriate to choose a new coordinate system 
rotating the original one used in (1) by the angle 9 in the xy plane (^(J. The rotation angle 
9 is determined by the requirement < S x ' >= in the new coordinate system. 
The retarded Green's function to be calculated is defined as 

Gy(t)=<S+(t);Sr(0)», (2) 

where S + and S~ are the spin-| operators in the rotated system. On introducing the two- 
dimensional Fourier transform G ninj (k\\,u}), one has the following form: 

< St; SJ^ > w = ^exp(zk||(rj - r^Gn^ik^uj), (3) 
k n 

where N' is the number of sites in any of the lattice planes, and rij represent the position 
vectors of site i and the layer index, respectively, kn = (k x ,k y ) is a two-dimensional wave 
vector parallel to the surface. The summation is taken over the Brillouin zone. 

As a result the equation of motion for the Green's function in Eq. (JSJ) of the ferroelectric 
thin film for T < T c has the following matrix form: 

H(w)G(k||,w) = R, (4) 
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where H can be expressed as: 

( u- V l 



with 



H 



o 



ki 

uo-V 2 k 2 
k 3 u-V 3 k 3 



\ 



k N u — Vn J 



h 

Vi 

V 2 
Vn 



J b a n sm 2 9 n , n = l,...,N, 

2Q S sin 9i + -a\ J s cos 2 Q\ — —7-^ sin 2 #17(^1) + J b a 2 cos 2 9 2 , 



1 

— < 

2 



2^5 sin 9 2 + 2 <T2 ^ b cos2 ^ 2 s * n2 ^"K^ll) + ^ sC71 cos2 0+ J&°3 cos2 ^3> 



= 2fi ra sin # n + -a n Jb cos 0. 



sin 2 6* n 7(k||) + J n _icr n _i cos 2 6'„_ 1 



V 
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iV 



2f2 s sin 9 N + -a N J s cos 2 9jy — ^ 



sin 2 0jv7(k|i) + Jb&N-i cos 2 0jv-i, 



7(ky) = -(cos(k x a) + cos(k y a)). 



Here we have introduced the notations J\ = Jtv = J s , Jn = Jb for n = 2, 3, 4, N — 1, 
fii = Qjy = Q s , Q n = Q b (n=2,3,4,...,N-l), J = Jn+i = 0. The quantity cr(T) is the relative 
polarization in the direction of the mean field and is equal to 2(S Z '). For the rotation angle 
9 we have the following two solution in the generalized Hartree-Fock approximation: 

7T 



1. cos9 = 0, i.e. 9 

2. sin 9 = — = — , 

a J a 



if T>T C - 
if T<T C . 



In order to obtain the solutions of the matrix equation (j%jL we introduce the two- 
dimensional column matrices, G m and R m , where the elements are given by (G n ) m = G mn 
and (R n )m = Cndmn, so that Eq. (JH) yields 



H(o;)G n = Rn 



(5) 



From Eq. (jSj), G nn (u) is obtained as: 



The quantity | H nn (u) | is the determinant made by replacing the n-th column of the 
determinant | H(lo) | by R n . The poles u n of the Green's function G nn (uj) can be calculated 
by solving | H{oj) |= 0. 

The relative polarization of the of the n-th layer is given by 



cr n Jn ^ 1-0-5 sin 2 fl n 7(k||) ^^y 1 (7) 
k n 

The last equation has to be calculated numerically. Due to the assumption of symmetrical 
surfaces, there are | N layer polarizations, which have to be solved self-consistent ly. In order 
to obtain the dependence of the Curie temperature Tc on the film thickness N, we let all 
a's be zero in Eq. ((7|) and solve the expressions self-consistently. 



III. NUMERICAL RESULTS AND DISCUSSION 



In this section we present the numerical results based on our theoretical calculations 
taking the typical the bulk parameter Jf, = 495K and Qb = 20K. We have calculated the 
temperature dependence of the polarization defined by 

a=< S z >=^$>n(T) cos0 n , 

n 

and the spin-wave frequencies of thin films for different values of the exchange interaction 
constants. One has to solve self-consistently the \ N coupled Eqs. (J7J) to obtain the layer 
polarization from the surface to half of the film (the other half is symmetric due to the 
assumption of identical surfaces). The numerical results expose some interesting and novel 
characteristics in the polarization and the spin- wave energies in comparison to the case of FE 
thin films without defects. To characterize the complete ferroelectric system both quantities, 
the polarization a and the spin-wave energy u , respectively are averaged over the iV layers. 
The results for film thickness iV = 9 and different exchange interaction parameters in the 
defect layers Jd are presented in Fig. ^ and Fig. |21 Let us firstly consider the case that 
the middle layer is defect layer. For example, it can be originated by localized vacancies 
or impurities with smaller radii and larger distances between them in comparison to the 
host material. It is reasonable to assume that the exchange interaction Jd is smaller than 
the value of the bulk interaction = 495K and has the value J 5 = J d = 300K (compare 
Figs. H 121 curve 1). The polarization (respectively the spin- wave energy) is smaller than 
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the case without defects, J5 = Jb (see Figs. [TJ El curve 2). The polarization decreases with 
increasing temperature to vanish at the critical temperature Tc of the thin film. The critical 
temperature decreases due to the smaller J d value. Pontes et al. [3] carried out dielectric 
and Raman spectroscopy studies and obtained that with addition of Sr to PbTiC>3, the 
phase-transition temperature decreases with the increasing Sr concentration. 

For the case where J5 = Jd = 1000 A" (Figs. [TJ curve 3), i.e. Jd is larger than the value 
of the bulk interaction constant J&, (for example when the impurities have a larger radius 
compared with the constituent ions) the polarization (respectively the spin-wave energy) is 
larger than the case without defects, J5 = Jb- Tc of the film is enhanced in comparison to 
the bulk value without defects due to the presence of larger J d values. This is the opposite 
behavior compared to the case of Jd = 300 AT, J d < Jb- Of course, in the case of both 
strong and small Jd values, thin film polarizations tend to the 3D one when n increases, 
since the excitation spectrum varies from 2D to 3D behavior with thickness. The second 
case, where Jd > Jb could explain the experimentally obtained increase in polarization 
the substitution of impurities, such as Nd, V, W and Nb in the layers of BTO films Q, 
or by increasing Ca contents in SBT thin films 11]. 

From Figs. ^ 121 on e can observe that the polarization, the spin-wave energies and the 
critical temperature of the FE phase transition are increased or decreased due to different 
exchange interactions in the defect layers. Our results are in qualitatively agreement with 
the experimental data of Noguchi et al. Q]. They have studied SBT films (Tc = 295°C) by 
the substitution of rare earth cations of La, Ce, Pr, Nd and Sm as well as Bi at the A site (Sr 
site) with Sr vacancies and have shown, that La modification induces soft behavior (lower 
E c and lower a r ), while a large amount of Nd and Sm substitution results in a very high E c 
(hard), as a result of defect engineering of both Sr and oxide vacancies. For not only SBT, 
but also other bismuth layer-structured FE (BLSF's) T c is strongly influenced by of A-site 
cations, and BLSF's with smaller A-site cations (Ca 2+ ) tend to show a higher Tc (420°C). 
The same amount of larger Ba 2+ brings about a relaxation of FE distortions and resulted in 
a decrease in Tc to 120°C The substitution of La led to a marked decrease in Tc to 180°C 
(x=0.5), because the induced A-site vacancies weaken the coupling between neighboring 
BOq octahedral P, 3|). This result corresponds in our calculations to the case of smaller 
values of the interaction constant in the defect layer Jd <Jt- For La- modified PbTiC>3, Tc 
decreased significantly, too, with an increase in La content I2J]. For Bi-SBT Tc rose strongly 
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to 405°C (:r=0.2) j 1 51 ] - The increase in Tc by Bi substitution is an opposite tendency to 
that in the case for La-SBT. The dominant role plays here the bonding characteristics with 
oxide ions. The influence of the orbital hybridization on Tc is very large, and Bi substitution 
resulted in a higher Tc Q|. This experimental result can be described qualitatively good 
using in our model exchange interaction parameters for the defect layer Jd > Jb- 

In Figs. ID is shown the influence of the surface exchange interaction constant J s on the 
polarization and spin-wave energy by J4 = Jd = const for a thin film with N = 7 layers. 
It can be seen that the two quantities are enhanced with increasing of J s . The critical 
temperature Tc increases, too. 

The polarization is depending, too, on that how many inner layers are defect. This is 
shown in Figs. [31 [U] for a thin film with N = 7 layers. In Fig. [S] is demonstrated that for 
Jd = 200K with increasing of the number of defect layers (curve 1 - J4 = Jd = Jb, curve 2 
- J 4 = Jd = 200K, curve 3 - J 3 = J 4 = J 5 = Jd = 200K) the polarization and the critical 
temperature decreases. The opposite behavior can be seen in Fig. [HI where Jd = lOOOK. 
With increasing of the number of defect layers the polarization and Tc increase. The spin- 
wave energies are reduced in the first case and enhanced in the second case. 

We have also studied a ferroelectric thin film whit N = 9 layers where one layer, the 
middle one is defect and the spins are ordered antiparallel, the defect layer is antiferroelectric, 
J 5 = J d = -1000K (Fig. [7[ curve 2) and J 5 = J d = -300K (Fig. El curve 2). The 
polarization and Tc of the film are smaller in comparison with the case where the middle layer 
is defect, but with ferroelectric order. The spin- wave energy of the thin film is also reduced 
when J5 = Jd = — 1000-fT. Such a sublattice model with ferroelectric and antiferroelectric 
layers where one or more layers are defect will be studied in a next paper. 



IV. CONCLUSIONS 



In conclusions, based on a modified transverse Ising model and using a Green's function 
technique, the polarization, the spin-wave energies and the phase transition temperature for 
ferroelectric thin films with structural defects are calculated. The dependence on temper- 
ature, film thickness and interaction constants is discussed. It is shown that defect layers 
in FE thin films, layers with impurities or vacancies, with dislocations, can induce strong 
increasing or decreasing of the polarization, the spin-wave energies and of the critical tem- 
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perature of FE phase transition due to different exchange interactions in the defect layers. 
The results are in good agreement with the experimental data for different ferroelectrics. 
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T[K] 

FIG. 1: Temperature dependence of the polarization a for a FE thin film with = 495K, Qb = 
20K, J s = 900K, Q s = Q b , N=9 and different ^-values: (I) J d =300 K, T c = 201K (□), (2) 495K, 
T c = 208K (o), (3) 1000K, T c = 224K (o). 
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FIG. 2: Temperature dependence of the spin-wave energy w[cm _1 ] for a FE thin film with J b = 
495K, n b = 20K, J s = 900K, tt s = Q b , N=9 and different J d -values: (1) J d = 300i-C (□), (2) 
495^(0), (3) 1000K(o). 
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FIG. 3: Temperature dependence of the polarization a for J b = 495X, £l b = 20K, J d = 1000i"C, 
n s = n b , N=7 and different J s -values: (1) J s = 200K, T c = Y&K (□), (2) 500, T c = 184K (o), 
(3) 1000K, T c = 237 K (o). 
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FIG. 5: Temperature dependence of the polarization a for J b = 495-fT, f2& = 20K, J s = 900K, 
Q s = Q b , N=7 and different defect layers: (1) J 4 = J b = 495K, T c = 220K(n), (2) J 4 = J d = 
200K, T c = 20OK'(o), (3) J 3 = J 4 = J 5 = J d = 200K, T c = 50K(o). 
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FIG. 6: Temperature dependence of the polarization a for J b = 495K, Q b = 20K, J s = 900K, 
n s = n b , N=7 and different defect layers: (1) J 4 = J b = 495K, T c = 220(d), (2) J 4 = J d = 
1000K, T c = 229(o), (3) J 3 = J 4 = J 5 = 1000K, T c = 259K(o). 
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FIG. 7: Temperature dependence of the polarization a with J b = 495K, £l b = 20K, J s = 900K, 
Q s = Q b , N=9 and an antiferroelectric defect layer: (1) J d = W00K, T c = 223K(n), (2) 
-lOOOtf, T c = 168«"(o). 
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FIG. 8: Temperature dependence of the polarization a with = 495K, fib = 20K, J s = 900K, 
fl s = Q b , N=9 and an antiferroelectric defect layer: (1) J d = 300K, T c = 20L?C(n), (2) 
-300K, T c = 190A"(o). 
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